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^ ' NA60 is a fixed-target experiment at the CERN SPS which measured dimuon production in 



nucleus-nucleus and proton-nucleus collisions. In this paper we report on a precision mea- 
surement of low-mass muon pairs in 158 AGeV indium-indium collisions. A significant excess 
of pairs is observed above the yield expected from neutral meson decays. The excess can be 
, isolated by subtraction of expected sources, thanks to the excellent mass resolution and large 

\^ • sample size. 

o ■ 

Q \ The ultimate goal of heavy ion collisions is the detection of signatures of a phase transition 

from hadronic matter to a deconfined and/or chirally restored medium. One of the main ob- 
^ . servation from the heavy ion program at the CERN SPS is the e xcess of low mass dimuons in 

pH \ nucleus-nucleus collisions reported by the CERES experiment^^. This is usually understood in 

terms of in-medium modification of the p meson, which could convey information about chiral 
symmetry restoration: changes in mass and/or width of the p are predicted as the chiral phase 
t raiis ition is approached. Most of the models agree in predicting an increase in the width of the 
■ piniS while several different predictions were formulated for its mass. The CERES data could 

be described by several significantly different models, so that an unambiguous interpretation of 
this result was not possible. A clear discrimination between different theoretical models requires 
good statistics and mass resolution, together with a good signal to background ratio. These 
requirements are met by the NA60 experiment. 

NA60 is a fixed target experiment which measures dimuon production in nucleus-nucleus and 
proton- nucleus collisions at the CERN SPS. The apparatus is composed of 4 main detectors: (^ij a 
muon spectrometer, which tracks muons and provides the trigger signal to the experiment; (ii) a. 
Zero Degree Calorimeter, for centrality measurements; (in) a silicon vertex tracker, which tracks 
particles in the vertex region; (iv) a silicon beam tracker, which tracks incoming beam particles 
before they hit the target. A complete description of the apparatus can be found elsewhere!^. 
Here we only briefly describe the detector concept. Experiments measuring dimuons usually put 
a hadron absorber in front of the muon spectrometer, to make sure only muons can reach this 
detector. This is also the main limiting factor, because of energy loss fluctuations and multiple 
scattering. NA60 overcomes this limitations by means of a vertex trackerEl, placed before the 
hadron absorber. Muons reconstructed in the muon spectrometer are extrapolated back to the 
vertex region and matched to vertex tracks, using both angles and momenta. This leads to 
a much improved vertexing (the origin of the muons can be accurately determined) and mass 
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Figure 1: Mass spectrum after combinatorial background (dashed) and fake matches (dot-dashed) subtraction 
(left). Isolation of excess by subtraction of expected sources: total data (open symbols), expected sources (contious 
line), difference data (thick triangles), sum of cocktail sources and difference data (dashed). 



resolution (from roughly 80 MeV to 20 MeV at the (j) mass). 

NA60 physics programme includes topics from the three mass region in which the dimuon 
invariant mass spectrum is customarily divided In this paper results on the Low Mass Region 
(LMR, M < 1.2 GeV) are reported. The data were collected in indium-indium collisions in 
2003; the sample analyzed in this work consists of roughly half of the 230 million events put 
on tape. The data were selected requiring only one vertex in the target region, to avoid re- 
interactions and pile-up events. The tracks reconstructed in the muon spectrometer were then 
matched to vertex tracks as stated above, allowing to extract the invariant mass spectrum of 
the resulting muon pairs. This is affected by two main sources of background: the combinatorial 
background and the fake matches, which have to be subtracted. The combinatorial background 
is the contribution to the dimuon spectrum of uncorrelated pairs coming from decays of vr and 
K. This is estimated with an event mixing technique: two muons from different events are 
randomly paired to build an invariant mass spectrum, which is uncorrelated by construction. 
This technique automatically takes into account all experimental details, provided the muons 
are taken with the correct normalization. The NA60 apparatus not only triggers opposite sign 
pairs (fi^fi^) but also like sign pairs {fJ,~ and which are made of uncorrelated muons 

only. The real and mixed spectra for the like-sign should then be identical. This can be used 
to asses the quality of the estimated background. The accuracy was found to be of the order of 
1% over several orders of magnitude^. As for fake matches, they are mistakes produced by the 
matching algorithm: a muon track can be matched to the wrong track in the vertex telescope, 
giving rise to a distorted spectrum. These can be estimated in two different ways. The first 
approach is an overlay Monte Carlo technique, where a Monte Carlo dimuon is reconstructed on 
top of a real event, allowing to check the probability of getting a fake match. The second method 
is an event mixing method, which extracts the probability for fake matches from data alone. 
The basic idea is to match the tracks in the muon spectrometer from one event to the vertex 
tracks of a different event. All the matches obtained in this way are fake by construction. This 
technique is more complicated but more rigorous. It is not used in the present analysis, though 
is crucial for other NA60 physics 

topicJSl. 

The two methods agree within 5%. Fake matches yield 
is less than 10% of the combinatorial background. The quality of the data after combinatorial 
background and fake matches subtraction is very good (Fig.0 left): the signal/background ratio 
is 1/7 and the mass resolution at the (p peak is 23 MeV. The analysis of the LMR data was 
done in 4 centrality bins, selected using charged particle multiplicity density: peripheral (4-30), 
semiperipheral (30- 110), semicentral (110-170) and central (170-240). The signal to background 
ratio, from peripheral to central is: 2, 1/3, 1/8, 1/11. The acceptance of NA60 extends all the 
way down low masses and pT, as opposed to previous dimuon experiments. This is important 
because the excess observed by CERES in the dielectron channel is stronger at low pt. 
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Figure 2: Excess mass spectra. The cocktail (solid) and the level of uncorrelated charm decays (dashed) are 

shown for comparison. The errors are purely statistical. 



The first step in the analysis was the study of the peripheral bin. In order to extract particle 
ratios in full phase space, data were fitted with expected sources: r] — > fifi, to — > /i/i, (p — > fifi, 
p — > /u/x, u — > HUTT^, T] — > Hfi'y, T]' — > nfi'y. The fit was done in 3 pT bins. The quality of the 
fit is very good. The results are found to be nearly independent of and in agreement with 
expectations, confirming that the acceptance is well under control. The particle ratios obtained 
from the fits without selection, corrected for acceptance and extrapolated to full phase space 
are rj/co = 0.88 ±0.03, (p/oj = 0.094 ±0.004 and p/oj = 1.62 ±0.10. The errors quoted are purely 
statistical. The value for ij/uj agrees with previously published results^. The value for (f>/uj is 
higher than measured in pp and pA collisions, as expected from strangeness enhancement. The 
value for p/u) is higher than expected from pp or pA, probably due to the fact that even in 
peripheral indium-indium collisions some hadronic medium is created, where the p is produced 
by vrvr annihilation. 

The more central bins could not be studied by a simple fitting procedure, because of the 
unknown source reported by CERES. Data were compared to a conservative cocktail, defined 
as follows: (i) the yields of (j) and to is fixed from the 2-body-decay peaks under the assumption 
that the underlying continuum is smooth (i.e. so as not to create any bump or dip in the 
spectrum after subtraction of these sources); (ii) the yield of the t] is fixed from the Dalitz 
channel saturating the very low mass part of the spectrum; (in) the yield of the p is fixed by 
assuming a p/u; ratio of 1.2 (as expected from previous pp and pBe measurements and from 
statistical models). The data show a clear excess when compared to this cocktail, becoming 
more pronounced with increasing centrality. Thanks to the good mass resolution, the excess 
could be isolated by subtraction of the known sources!^ (Fig-E right), namely the conservative 
cocktail described above, with the exception of the p: since the modification in the spectrum 
is generally understood as due to in-medium modification of this particle, it does not make 
sense to subtract from the spectrum its vacuum line-shape. The difference spectrum obtained 
in this way is a lower limit to the excess, as a consequence of the procedure adopted to fix yields 
(for example, there is no excess at very low masses by construction). Fig. [21 shows the excess 
as a function of centrality. For comparison, the cocktail p and DD are also plotted. A clear 
excess rising with centrality and centered at the nominal p pole is visible. The excess spectrum 
is robust against systematical uncertainties: both the combinatorial background and the fake 
matches are flat as a function of mass. A 1% error on combinatorial background or a 5% error 
on fake matches would not change the overall features of the spectrum. Changing the procedure 
to fix the r] normalization, for example by requiring that it accounts for 90% of the data at very 
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Figure 3: Comparison of the semi-peripheral excess mass spectrum to model predictions. The open data points 
show the difference spectrum obtained decreasing the -q yield by 10%, reflecting the systematic uncertainty in 
this normalization. Left: Cocktail p (thin solid), unmodified p (dashed), in-medium broadening (thick solid), 
dropping mass * (dashed-dotted). Right: Comparison to models of Ref. ^'^ (solid) and Ref. ^ (dashed-dotted). 



low masses, has a small effect on the lowest part of the spectrum only, while the yield above 
0.4 GeV is unaffected (Fig.©. 

The excess was compared to theoretical models available before the NA60 datJ3^. Dropping 
mass scenarios seem to be ruled out, while broadening models are strongly favored as can be seen 
in the left part of Fig. |21 (theoretical models are normalized to data in the range M < 0.9 GeV). 
None of the existing models was able to describe the data in full detail. The qualitative features 
of the spectrum are consistent with p production via tttt annihilation in the framework of a 
hadron many-body theory. The prediction, however, fails to describe the excess spectrum for 
M > 0.9 GeV. An updated version of this mode^^, implementing an improved fireball dynamics 
and 47r processes (including a contribution from vector-axialvector mixing) is able to describe 
the spectrum, even in absolute term (Fig. fright, solid curve). A different model, based on 
a similar idea but with a different technical implementation, considering on lv a hot pion gas 
without baryons (which are an essential ingredient in the previous approach) is also able to 
describe the data in absolute terms (Fig. 01 right, dashed-dotted curve). The An processes are not 
included in this model. Both models require a continuum contribution from a partonic spectral 
function, which is much stronger in the second one and fills the high mass region. A conclusive 
discrimination between different scenarios may come from the study of the dependence of the 
excess data. Attempts to theoretically understand NA60 data are progressing steadily and more 
approaches, different from the ones briefiy outlined here, are appearing in the literature'^. 
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